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• Abstract 

(N 

We investigate the in-medium modification of the charmonium break-up process 
due to the Mott effect for light (n) and open-charm (D, D*) mesons at the chi- 
ral/deconfinement phase transition. A model calculation for the process J/ip + ir — > 
D + D* + h.c. is presented which demonstrates that the Mott effect for the D-mesons 
■ leads to a threshold effect in the thermal averaged break-up cross section. This effect 

is suggested as an explanation of the phenomenon of anomalous J/ ip suppression in 
the CERN NA50 experiment. 
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1 Introduction 

S3' 



Recent results of the CERN NA50 collaboration on anomalous J/ip suppres- 
sion [1] in ultrarelativistic Pb-Pb collisions at 158 AGeV have renewed the 
quest for an explanation of the processes which may cause the rather sud- 
den drop of the J/ip production cross section for transverse energies above 
E T ~ 40 GeV in this experiment. An effect like this was predicted as a signal 
for quark gluon plasma formation [2] due to screening of the cc interaction. 
Soon after that it became clear that for temperatures and densities just above 
the deconfinement transition the Mott effect for the J/ip does not occur and 
that a kinetic process is required to dissolve the J / ip [3] in a break-up process 
by impact of thermal photons [4], quarks [5], gluons [6] or mesons [7,8]. 

In this paper, we suggest that at the chiral/deconfinement phase transition 
the charmonium break-up reaction cross sections are critically enhanced since 
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the light and open-charm mesonic states of the dissociation processes become 
unbound (Mott effect) so that the reaction thresholds are effectively lowered. 
We present a model calculation for the particular process J /if) + n — > D + 
D* + h.c. in a hot gas of resonant (unbound but correlated) quark- ant iquark 
states in order to demonstrate that the Mott dissociation of the final states 
(D-mesons) at the chiral phase transition leads to a threshold effect for the 
in-medium J /if) break-up cross section and thus the survival probability. 



2 In-medium modification of charmonium break-up cross sections 

The inverse lifetime of a charmonium state in a hot and dense many-particle 
system is given by the imaginary part of its selfenergy 

r- 1 (p)=r(p) = S>(p)-S < (p). (1) 




Fig. 1. Diagrammatic representation of the complex selfenergy for the J /if) due to 
break-up in (off-shell) D, D* pairs by impact of (off-shell) pions from a hot medium. 

In the Born collision approximation for the dominant process in a hot gas of 
pion-like correlations, as shown in Fig. 1, we have [9] 

= / // (27r) 4 <W;^ 2 \M\ 2 G>(P') G< 2 (p 2 ) , (2) 

p> Pi Pi 

where the thermal Green functions Gf(p) = [1 + fi(p)]Ai(p) and Gf{j>) = 
fi(p)Ai(p) are defined by the spectral function Ai{p) and the distribution 
function fi(p) of the bosonic state i; with the notation J p = J (f^r, ^p,p';pi,p 2 = 
S(p + p' - pi -p 2 )- 
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In the low density approximation for the final states (/^(p) ~ 0), one can 
safely neglect S < (p) so that 

T ~ 1{p) = 1 1 J(^)%p'-, P up,\M\ 2 U(P') MP') A Dl ( Pl ) A D2 ( P2 ). (3) 

p' Pi P2 

With the differential cross section 

da = 1 \M(s,t)\ 2 

dt 16tt A(s, M2, s >) ' 1 J 

using s = (p + p') 2 > t = (p — pi) 2 , s' = p' 2 and A(s, M 2 , s') = [s — (M^ + 
v/s 7 ) 2 ] [s - (M^ - Vs 1 ) 2 } = 4 w 2 cl [p 2 + M|] [p /2 + s'] one can show that the J/ip 
relaxation time in a hot pion as well as pionic resonance gas is given by 




where depending on the properties of the medium the pion spectral function 
describes either qq bound states or resonant (off-shell) correlations. The in- 
medium break-up cross section is given by 

cr*(s) = J dsi ds 2 A Dl (si) A D2 (s 2 ) a(s;s 1 ,s 2 ) ■ (6) 

Note that there are two kinds of medium effects due to (i) the spectral func- 
tions of the final states and (ii) the explicit medium dependence of the matrix 
element M.. In the following model calculation we will use the approximation 
a(s; si, s 2 ) ~ o" vac (s; s±, s 2 ) justified by the locality of the transition matrix M. 
which makes it rather inert against medium influence. 



3 Model calculation 

The quark exchange processes in meson-meson scattering can be calculated 
within the diagrammatic approach of Barnes and Swanson [10] which allows 
a generalization to finite temperatures in the thermodynamic Green function 
technique [11]. This technique has been applied to the calculation of J/ip 
break-up cross sections by pion impact in [8]. The approach has been ex- 
tended to excited charmonia states and consideration of rho-meson impact 
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recently [12]. The generic form of the resulting cross section (given a band of 
uncertainty) can be fit to the form 

a™(s; M 2 Di ,M 2 D2 ) = a ln(s/s ) exp(-s/A 2 ) , s > s , (7) 

where s = (M Dl +M D2 ) 2 is the threshold for the process to occur, a = 7.5-10 9 
mb and A = 0.9 GeV. 

Recently, the charmonium dissociation processes have been calculated also in 
an effective Lagrangian approach [13,14], but the freedom of choice for the 
formfactors of meson-meson vertices makes predictions uncertain. The devel- 
opment of a unifying approach on the basis of a relativistic confining quark 
model is in progress [15] and will remove this uncertainty by providing a deriva- 
tion of the appropriate formfactors from the underlying quark substructure. 

The major modification of the charmonium break-up process which we expect 
at finite temperatures in a hot medium of strongly correlated quark- ant iquark 
states comes from the Mott effect for the light as well as the open-charm 
mesons. At finite temperatures when the chiral symmetry in the light quark 
sector is restored, the continuum threshold for light-heavy quark pairs drops 
below the mass of the D-mesons so that they are no longer bound states 
constrained to their mass shell, but become rather broad resonant correlations 
in the continuum. This Mott effect has been discussed within relativistic quark 
models [16] for the light meson sector but can also be generalized to the case 
of heavy mesons [17]. Applying a confining quark model [18] we have obtained 
the critical temperatures T^ ott = 110 MeV, T^ ott = 140 MeV and T^ ott = 150 
MeV [19]. 

In order to study the implications of the pion and D-meson Mott effect for 
the charmonium break-up we adopt here a Breit-Wigner form for the spectral 
functions 



l{ } vr (s - Mf(T)y + T 2 (T)M 2 (T) ' { } 

which in the limit of vanishing width Tj(T) — > goes over into the delta 
function 5(s — M 2 ) for a bound state in the channel i. The width of the pions 
as well as the D-mesons shall be modeled by a microscopic approach. For our 
exploratory calculation, we adopt here 

r W|D (T) = c (T - Tffi ) 0(T - T^) , (9) 



where the coefficient c = 2.67 is assumed to be universal for the pions and 
D-mesons and it is obtained from a fit to the pion width above the pion Mott 
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temperature, see [20]. For the meson masses we have M^^iT) = M n> D + 
0.75 T Wt £)(T). The result for the in-medium J /if) break-up cross section (6) is 
shown in Fig. 2. 




Fig. 2. Energy- and temperature dependent in-medium J /if) break-up cross section 
for pion impact. Thresholds occur at the Mott temperatures for the open-charm 
mesons: T^, ott = 110 MeV, T^ ott = 140 MeV. 

With M D , = 2.01 GeV and M D = 1.87 GeV follows for the threshold s = 
15.05 GeV 2 . At a temperature T = 140 MeV, where the D -meson can still 
be considered as a true bound state, the _D*-meson has already entered the 
continuum and is a resonance with a half width of about 80 MeV. Due to 
the Mott effect for the open-charm mesons (final states), the charmonium 
dissociation processes become "subthreshold" ones and their cross sections 
which are peaked at threshold rise and spread to lower onset with cms energy. 
This is expected to enhance strongly the rate for the charmonium dissociation 
processes in a hot resonance gas. 



4 J /if) dissociation in a hot "pion" gas 

We calculate the inverse relaxation time for a J /if) at rest in a hot gas of pions 
(below T^ ott ) and pion-like qq correlations (above T^ ott ) by specifying Eq. 
(5) for this case 
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r ^ T )= fw^sf ds^MW^T) - , |p/| v*(s) (10) 

= <a*v Tcl >n T (T) , (11) 

with the dispersion relation E^(jp',s n ) = J p' 2 + s n , the thermal Bose distri- 
bution function f n (p',s n ;T) = 3 {exp[i? 7r (p / , s n )/T] — 1} _1 and the particle 
density n n (T) for the "pions". The cms energy of the "pion" impact on a J/if> 
at rest is s(p'; s^) = s„- + M 2 + 2M^,E 7r (p / , s^). 

The result for the temperature dependence of the thermal averaged J j ip break- 
up cross section < a*v rc \ > is shown in Fig. 3. This quantity has to be com- 
pared to the nuclear absorption cross section for the J / ip of about 3 mb which 
has been extracted from charmonium suppression data in p-A collisions [21]. 




T [MeV] 

Fig. 3. Temperature dependence of the thermal averaged in-medium J ftp break-up 
cross section for different charm quark masses. 



It is remarkable that it is practically negligible below the D-meson Mott tem- 
perature T]^ ott = 110 MeV but comparable to the nuclear absorption cross 
section above the chiral/deconfmement temperature of T cr it ~ 150 MeV. It 
is obvious that the transition from D-meson bound states to unbound light- 
heavy quark correlations is responsible for the strong increase by one to two 
orders of magnitude. Note that in this calculation the Mott effect for the pion 
(initial state) above T™ ott has been included, but does not alter the result 
obtained previously [19] in a calculation neglecting this effect. 

Therefore we expect the in-medium enhanced charmonium dissociation pro- 
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cess to be sufficiently effective to destroy the charmonium state on its way 
through the hot fireball of the heavy-ion collision and to provide an explana- 
tion of the observed anomalous J/ip suppression phenomenon [1]. A detailed 
comparison with the recent data from the NA50 collaboration requires a model 
for the heavy-ion collision. The effective in-medium break-up cross section for 
the J/ifi derived in this work provides an input for all calculations which use 
this quantity, e.g. Glauber-type models [19,22-26], more detailed calculations 
based on a parton cascade model [27] or molecular dynamics [28]. 



5 Summary and Outlook 

In this letter we have presented an approach to charmonium break-up in a hot 
and dense medium which is applicable in the vicinity of the chiral/deconfine- 
ment phase transition where mesonic bound states get dissolved in a Mott- 
type transition and should be described as resonant correlations in the quark 
plasma. This description can be achieved using the concept of the spectral 
function which can be obtained from relativistic quark models in a system- 
atic way. The result of an exploratory calculation employing a temperature- 
dependent Breit-Wigner spectral function for light and open-charm mesons 
presented in this paper has demonstrated that heavy-flavor dissociation pro- 
cesses are critically enhanced at the QCD phase transition and could represent 
the physical mechanism behind the phenomenon of anomalous J/ip suppres- 
sion. 

In subsequent work we will relax systematically approximations which have 
been made in the present paper and improve inputs which have been used. 
In particular, we will investigate the off-shell behaviour of the charmonium 
break-up cross section in the vacuum (7) and calculate the spectral functions 
(8) at finite temperature within a relativistic quark model. Dyson- Schwinger 
equations provide a nonperturbative, field-theoretical approach which has re- 
cently been applied also to heavy-meson observables [29] and have proven 
successful for finite-temperature generalization [30,31]. Further intermediate 
open-charm states can be considered; the states in the dense environment 
should include rho mesons and nucleons besides of the pions which all can be 
treated as off-shell quark correlations at the QCD phase transition. 

In future experiments at LHC the charm distribution in the created fireball 
may be not negligible so that the approximation foXp) ~ has to be re- 
laxed. In this case, one has to include the gain process (DD annihilation) 
encoded in the function. In comparision to previous investigations [32,33] 
the present quantum kinetic treatment contains Bose enhancement factors in 
the Gf functions which modify the charm equilibration process. 
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